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Abstract 

The spectral structure of the Laplacian-Beltrami operator (LBO) on manifolds has been widely used 
in many applications, include spectral clustering, dimensionality reduction, mesh smoothing, compression 
and editing, shape segmentation, matching and parameterization, and so on. Typically, the underlying 
Riemannian manifold is unknown and often given by a set of sample points. The spectral structure of 
the LBO is estimated from some discrete Laplace operator constructed from this set of sample points. 
In our previous papers |13lI17| . we proposed the point integral method to discretize the LBO from point 
clouds, which is also capable to solve the eigenproblem. Then one fundmental issue is the convergence 
of the eigensystem of the discrete Laplacian to that of the LBO. In this paper, for compact manifolds 
isometrically embedded in Euclidean spaces possibly with boundary, we show that the eigenvalues and the 
eigenvectors obtained by the point integral method converges to the eigenvalues and the eigenfunctions 
of the LBO with the Neumann boundary, and in addition, we give an estimate of the convergence rate. 
This result provides a solid mathematical foundation for the point integral method in the computation 
of Laplacian spectra from point clouds. 
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1 Introduction 

The Laplace-Beltrami operator (LBO) is a fundamental object associated to Riemannian manifolds, which 
encodes all intrinsic geometry of the manifolds and has many desirable properties. It is also related to 
diffusion and heat equation on the manifold, and is connected to a large body of classical mathematics (see, 
e.g., [H]). In recent years, the Laplace-Beltrami operator has attracted much attention in many applied 
fields. The spectral structure of the Laplacian-Beltrami operator on manifold has been widely used in many 
applications, include spectral clustering, dimension reduction, mesh smoothing, compression, editing, shape 
segmentation, matching, parameterization, and so on [siiniiisiiii]. 

In this paper, we consider the following eigenproblem on a smooth, compact A:-dimensional submanifold 
M with the Neumann boundary condition 


A^u(x) = Xu{x), X £ VI 

f^(x) = 0, KGdM 


( 1 . 1 ) 


where is the Laplace-Beltrami operator on A4. Let g be the Riemannian metric tensor of M, which 
is assumed to be inherited from the ambient space that is, Ai isometrically embedded in R'^ with the 

standard Euclidean metric. If Ai is an open set in then A^ becomes standard Laplace operator, i.e., 

A _ 
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In [13], we have proposed a numerical method called the point integral method (PIM) to solve the above 
eigenproblem on manifolds. In the PIM, we only need a point cloud to discretize the manifold A4. In 
particular, we are given a set of points P = (pi, • • • , p„) sampling A4. In addition, we are given a vector 
V = (Vi, • • • , !/„)* where Vi is an volume weight of Pi in AI, so that for any Lipschitz function f on A4, 

f{:x.)dfi:xL can be approximated by /(Pi)^i- Here d/ix is the volume form of Ai. 

We consider the following discrete Laplace operator Lt^h 

Lt,hu{pi) = “ u{pj))Vj. (1.2) 

PjeP ^ ^ 

where R : M'*' —>■ ffi.'*' is a given kernel function and Ct = is a normalized constant. If set V) = ^, the 

discrete Laplace operator Lt^h becomes the well-known weighted graph Laplacian |5]. Denote R the primitive 
function of —R, i.e. -£:R{r) = —R{r). Now we consider the following generalized eigenproblem of Lt^h 

- Lt.,zr(p,) =Xj2CtR u(p^)V, (1.3) 

The purpose of the paper is to show the generalized eigenproblem (11.311 converges to the eigenproblem (11.111 
at a rate (t^^^ -I- jkjk+a) for the eigenvalues and (t^/^ -b ^k/l+ 2 ) for the eigenfunctions, provided the input 
data (H, V) is an /i-integral approximation of Ai (see Sectionfor its definition). Note the rate reported 
in this paper depends on the dimension k of Ai and may not be optimal. 

Following [nun], we bridge the LBO Am and the discrete Lapalce operator Lt^h using the following 
integral Laplace operator 

LtM(x) = ^ J^R (m(x) - M(y))d/iy. (1.4) 

We consider the solution operators associated with Lt and Lt^h, and show all solution operators are in 
fact compact and additionally the approximation errors of the solution operators in certain operator norms 
are bounded (see Theorem l4.3l and l4.4ll . This enables us to apply the results from spectral theory to obtain a 
convergence rate. Note that it is critical and also common in the numerical analysis to consider the solution 
operators [6] , instead of the Laplacians themselves which are not even bounded. Comparing to the result in 
m we improve the error estimations for both the trancation error and the stability of Lt- In particular, 
in both Theorem O and 15.31 we bound the approximation errors using the Sobolev norms, instead of the 
infinite norm as in HZ]. 

1.1 Related work 

The finite element method is popular in solving PDFs on manifolds. Dziuk m showed the FEM converges 
quadratically in norm and linearly in norm in solving the Poisson equations. The eigensystem of the 
LBO computed by the FEM converges linearly [Hiiniiin]- The FEM requires a mesh tesselating the domain 
and its performace depends heavily on the quality of meshes. However, it is not an easy task to generate a 
mesh of good quality, especially for an irregular domain [8], and becomes much more difficult for a curved 
submanifold [7]. 

We see that the discrete Laplace operator Lt^h becomes the weighted graph Laplacian when the volume 
weight vector is constant. In the presence of no boundary, Belkin and Niyogi [3] showed that the spectra of 
the weighted graph Laplacian converges to the spectra of . When there is a boundary, it was observed 
in [HIS] that the integral Laplace operator Lt is dominated by the hrst order derivative and thus fails to 
be true Laplacian near the boundary. Recently, Singer and Wu [TS] showed the spectral convergence in 
the presence of the Neumann boundary. In the previous approaches, the convergence analysis is based on 
the connection between the weighted graph Laplacian and the heat operator, and thus it is essential to use 
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the Gaussian kernel in those approaches. The convergence analysis in this paper is very different from the 
previous ones. We consider this convergence problem from the point of view of solving the Poisson equation 
on submanifolds, which opens up many tools in the numerical analysis for studying the graph Laplacian. 

The rest of the paper is organized as follows. In Section [2l we describe the basic assumptions and define 
the solutions operators we are working on. The main results are stated in Section [3l In Section |4l the proofs 
of the main results are given. The proofs depends on two theorems: one states that Tt converges to T in 
norm and the other states that Tt^h converges to Tt in norm, whose proofs are given in Section [5] and | 6 ] 
respectively. Finally, in Section [51 we conclude and point out a few direction for future research. 


2 Assumptions and Notations 

We follow the assumptions in m and use basically the same notations. 


2.1 Assumptions 

First we assume the function R : R+ —>■ R+ is and satisfies the following conditions: 

(i) R{r) = 0 for Vr > 1. 

(ii) There exists a constant Sq so that R{r) > 6q for Vr < 1. 

Second, we assume both Ai and dA4 are compact and C°° smooth. Consequently, it is well known that 
both Ai and dA4 have positive reaches. 

Finally, we assume the input data (T*, V) is h-integral approximation of Ai, i.e.. 

For any function / G C^(A4), there is a constant C independent of h and / so that 


'M 


2 = 1 


< C'/i|supp(/)|||/||ci, and 


( 2 . 1 ) 


where ||/||ci = ll/lloo + ||V/||oo and A| denotes the volume of X. 

Remark 2.1. The quadrature rule we considered in this paper is in a deterministic way. We assume that 
the sample points P converge to the underlying manifold At in the sense that > 0. Actually, the methods 
used in this paper also apply on random samples and the results will be reported in the near future. 


2.2 Notations 

To investigate the convergence of the problem 1 1.31 to the problem 1 1.1 1 of the Neumann boundary, we consider 
the solution operators T, Tt and Tt^h as follows. Denote the operator T : Lf{Al) — >■ H‘^{Ai) to be the solution 
operator of the following problem 

J Xmu{A) = /(x), X G 7W, 

\ |H(x)= 0 , xG^M. 

where n is the out normal vector of Ai. Namely, for any / G Lf{Ai), u = T{f) with u = 0 solves the 
problem ( 12 . 21 ) . 

Denote Tt : Lf{Ai) —>■ Lf{Ai) to be the solution operator of the following problem 

- 7 / ■Rt(x,y)(M(x) - u(y))dy = / i? 4 (x, y)/(y)dy. (2.3) 

t JM JM 

where = CtR and .Ri(x,y) = CtR Rir) = f.^°° R{s)ds. Namely, u = Tt{f) with 

Im u = 0 solves the problem ( 12 .3|) : 
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The last solution operator is Tt^t '■ C(A4) —>■ C(M) which defined as follows. 



(2.4) 


n 


where Wt,h{^) = Rt{^, and u = (rti, • • • , UnY with UiVi = 0 solves the following linear system 




One direct consequence is that T, Tt and Tt^h have following properties. 

Proposition 2.1. For any t > 0, h > 0, 

1. T,Tt are compact operators on H^{A4) into H^(Ai); Tt,Tt^h o,re compact operators on C^(M) into 


C\M). 


2. All eigenvalues of T,Tt,Tt^h real numbers. All generalized eigenvectors of T,Tt,Tt^h orre eigenvec¬ 
tors. 

Proof. The proof of (1) is straightforward. First, it is well known that T is compact operator. Tt^t is actually 
finite dimensional operator, so it is also compact. To show the compactness of Tt, we need the following 
formula. 



( 2 . 6 ) 


Using the assumption that R € , direct calculation would gives that that Ttu G . This would imply 


the compactness of Tt both in and C^. 

For the operator T, the conclusion (2) is well known. The proof of Tt and Tt^h are very similar, so here 
we only present the proof for Tt. 

Let A be an eigenvalue of Tt and u is corresponding eigenfunction, then 


LtTtU = XLtU 


which implies that 


^ ^ /m /m -^t(x,y)M*(x)M(y)dxdy 
XMW*W(Ltn)(x)dx 


where u* is the complex conjugate of u. 

Using the symmetry of Lt and .R(x,y), it is easy to show that A £ R. 

Let M be a generalized eigenfunction of Tt with multiplicity m > 1 associate with eigenvalue A. Let 


V = {Tt — A)™ ^u, w = {Tt — A)™ then v is an eigenfunction of Tt and 


Ttv = An, {Tt — X)w = v 


By applying Lt on both sides of above two equations, we have 
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Using above two equations and the fact that Lt is symmetric, we get 


0 = (w,XLtV- f i?t(x,y)u(y)dy\ 

\ JM / M 

= iXLtW- / Rti:x.,y)w{y)dy,v) 

\ Jm / M 

= {Ltv,v)^ > C\\v\\l 

which implies that {Tt — = v = 0. This proves that u is a generalized eigenfunction of Tt with 

multiplicity m — 1. Repeating above argument, we can show that u is actually an eigenfunction of Tj. □ 

Using the definition of T^Tt and Tt^h, it is easy to show that the eigen problems Tu = Xu, Tt^h{u) = Xu 
is equivalent to the eigen problems (HU) and (HU respectively. Namely their eigenvalues are reciprocal to 
each other and they share the same eigenspaces. 

Proposition 2.2. Let 0(u) denote the restriction ofu to the sample points P, i.e., Oiu) = {u(pi), • • • , m(p„))*. 


1. If a function u is an eigenfunction ofTt^h with the eigenvalue X, then the vector 0(u) is an eigenvector 
of the eigenproblem HU with eigenvalue 1/A. 

2. If a vector u is an eigenvector of the eigenproblem (11.31) with the eigenvalue A 7 ^ 0, then /a(u) is an 
eigenfunction of Tt^h with eigenvalue 1/A, where 


^a(u)(x) 


Spj eP Pi Spj GP Pj ^ 

P3)Vj 


3. 


A function u is the eigenfunction of the eigenproblem HU with the eigenvalue X 0 if and only if the 
function u is an eigenfunction of T with the eigenvalue 1/A. 


Using the above proposition, we only need to prove the eigenvalues and the eigenfunctions of Tt^h converge 
to the eigenvalues and the eigenfunctions of T. 


3 Main Results 

Let X be a complex Banach space and L : X —>■ X be a compact linear operator. The resolvent set p{L) is 
given by the complex numbers z G C such that z — L is bijective. The spectrum of L is cr{L) = C\p{L). It 
is well known that cr(L) is a countable set with no limit points other than zero. All non-zero value s in a{L) 
are eigenvalues. If A is a nonzero eigenvalue of L, the ascent multiplicity a of A — L is the smallest integer 
such that ker(A — L)°‘ = ker(A — L)“+^. 

Given a closed smooth curve T C p{L) which encloses the eigenvalue A and no other elements of cr(L), 
the Riesz spectral projection associated with A is defined by 

E{X,L) = ^J^iz-L)-^dz, ( 3 . 1 ) 

where i = is the unit imaginary. The definition does not depend on the chosen of T. It is well known 
that E{X,L) : X —X has following properties: 

1. E{X, L) o E{X, L) = E{X, L), L o E{X, L) = E{X, L) o L, E{X, L) o E{pi, L) = 0, if A 7 ^ p. 

2. E{X, L)X = ker(A — L)“, where a is the ascent multiplicity of A — L. 
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3. If r C p{L) encloses more eigenvalues Ai, • • • , A„, then 

- ■ ■ ,Xn,L)X = (B'i—i ker(Ai — L)“* 
where ai is the ascent multiplicity of Ai — L. 

The properties (2) and (3) are of fundamental importance for the study of eigenvector approximation. 

It is well-known that the eigenvalues of both T can be sorted as 

Ai < A 2 < • • • , < A„ < • • • < 0, 

where the same eigenvalue are repeated accroding to its multiplicity. 

Now, we are ready to state the main theorems. 

Theorem 3.1. A ssume the submanifold M and dM are C°° smooth and the input data (T*, V) is an h- 
integral approximation of Ai, Let Ai be the ith smallest eigenvalue ofT counting multiplicity, and A*’^ be the 
ith smallest eigenvalue ofTt^h counting multiplicity, then there exists a constant C such that 

|A« - A.I < C . 

and there exist another constant C such that, for any 4> € E(Xi,T)X and X = H^(Ai), 

u - < C . 

where = {A*’^ e o-{Tt^h) ■ j G h} and A = {j G N : Xj = Ai}. 

Using the above theorem, we know that the eigenvalues of Tt^h converge to the eigenvalues of T. From 
ProDOsition l2.21 the inverses of the eigenvalues of Tt^h and T are the eigenvalues of (jl.3p and (11.11) respectively. 
Then, the eigenvalues of (11.31) converge to the eigenvalues of dnj. The above theorem also states the 
convergence of the generalized eigenfunctions, namely the generalized eigenfunctions of /j to the generalized 
eigenfunctions of T. Since all generalized eigenfunctions of Tt h and T are indeed the eigenfunctions as proved 
in Propositions l2.ll we obtain the convergence of the eigenfunctions of Tt^h to that of T. Then the convergence 
of the eigenfunctions of (HJl) to that of dm) can be obtained by using the facts that Tt^h and (11.31) (T and 
(ll.l[) l share the same eigenfunctions as proved in Proposition [521 


4 Structure of the Proof 


To prove Theorem 13.11 we need some results in perturbation theory for spectral projection. First, we need 
the following theorem to obtain the convergence rate of the eigenvalues. 


Theorem 4.1. Let {X, || • ||x) be an arbitrary Banach space. Let S and T he compact linear operators on 
X into X. Let z G p(T). Assume 


\\T-S\\x< 


1 


Then z G p{S) and {z — S) ^ has the bound 


\\{z-S)-^\\x<2\\{z-T)-^\U. 


(4.1) 


(4.2) 
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Proof. For any x G X 


||(^-^)x|U > \\{z-T)x\\x-\\{T-S)x\\x 

> (||(z-r)-i|U-||T-5||x)IN|x 

> l\\iz-T)-^x\\x\\x 

Using this inequality and the assumption that S is compact operator, we have 2 € p{S) and 

\\{z-S)-^\\x<2\\{z-T)-^x- 


(4.3) 


(4.4) 

□ 


We also need the following theorem (e.g. m) to get the convergence rate of the eigenfunctions. 

Theorem 4.2. Let {X, || • ||x) be an arbitrary Banach space. Let S and T be eompaet linear operators on 
X into X. Let zq G C, zq ^ 0 and let e > 0 be less than \zo\, denote the cireumference \z — Zo\ = e byV and 
assume F C p{T). Denote the interior of T by U. Let ax = U (1 a(T) ^0. as = U r\a{S). Let E{as,S) 
and E(aT,T) be the corresponding spectral projections of S for as and T for ax, i-c. 


Assume 


Then, we have 


E{as,S) = ^ ^(z - 5 )-Mz, E{ax,T) = ^ ^(z - T)-Mz. 


||(r- 5')5'||x < min- 
z^r 


(z-T)- 




(4.5) 


(4.6) 


(1) . Dimension E{as, S)X = E(ax,T)X, thereby as is nonempty and of the same multiplicity as ax- 

(2) . For every X G E{ax,T)X, 

||x - E{as, 5)x||x < — (||(T - 5)a;||x + ||x||x||(r - S)S\\x) - (4.7) 

Co 

where M = max^gr ||(-2 - T)~^\\x, co = min^gr ki¬ 


ln order to utilize the above two theorems from pertubation theory, we show the following three theorems 
which bound the approximation errors of the solution operators T,Tt,T^ h in certain operator norms. 

Theorem 4.3. Under the assumptions in \2.1[ there exist generic constants C > 0 and To > 0 only depend 
on j\4, such that 

\\T-Tt\\m<Ct^^^ (4.8) 


provides t <Tq. 

Theorem 4.4. Under the assumptions in \2.1l there exist generic constants C > 0 and Tq > 0 only depend 
on M., such that 

Ch 

im,.-T,||e.<^ (4.9) 

as long as t < Tq and <Tq. 
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Theorem 4.5. Under the assumptions in \2.1l there exist generic constants C > 0 and Tq > 0 only depend 
on j\4, such that 


WTtWrn < C, \\Tt,h\\oo < Ct-^^\ \\Tt,h\\c^ < (4.10) 

as long as t < Tq and <Tq. 

The above three theorems will be proved in Section O Section [5] and Section [7] respectively. 

Before we prove Theorem l3.ll we show the the estimates of \\{z — T)~^\\Hi-(M) \\{z — iii 

the following two lemmas, which are needed to invoke the results from pertubation theory. 

Lemma 4.1. Let T he the solution operator of the Neumann problem ()2.2p and z £ p{T), then 

||(z-r)-i||^i(;K) < max--, (4.11) 

^ ' nGN \Z — Xn\ 

where {A„}ngN the set of eigenvalues ofT. 


Proof. Suppose (fj , j £ N be the normalized eigenfunction of T corresponding to \j , j £ N. Then it is well 
known that is a orthonormal basis of 

For any x £ H^{A4), z £ p(T), first we can expand x over {<(>j}jGN to obtain 

OO 

^ = (4-12) 

1=1 


Then, we have 


iz-T)x\\H^ = 


'^Cj{z-T)(j)j 


1=1 


( OO 

= min| 2 -A„|||x|l//i 

kGN 


1/2 


J2cjiz- 

1=1 




/ OO 


1/2 


> min z — 
neN 




(4.13) 


□ 


Lemma 4.2. Let Tt he the solution operator of the integral equation (12.3p . For any z £ C\ i?(A„, cq) 
with rg > ||T —Ttll/ji, then 


IKz-Tt) ^||ci<max 


2|7W| 

|2,|^(fe+2)/4 


mm z — 

V riGN 


A„|- ||T-Tt||ffi 



Proof. For any x £ 


\\{z-Tt)x\\H^ > \\{z-T)x\\m-\\{T-Tt)x\\m 

ynGN J 


Then {z — Tt) ^ exists and 


Wiz-Tt) < ( min|z - A„| - ||r 



(4.14) 


(4.15) 


(4.16) 












For any u G C^{M), 


-1 


\\{z-Tt)-^u\\H^<[min\z-Xr.\-\\T-Tt\\H^) \M\\\u\\c^ (4.17) 

ynGN J 

where |A1| is the volume of the manifold Ai. 

On the other hand, let v = (z — Tt)~^u which means v = {u + Tfv)/z 

Ikllc^ < Adkllc^ + im^llcO 
1^1 

which proves that 

ll(^ - Tt)-^\\c^ < max (^mm \z - A„| - ||T - (4.18) 

□ 

Using the above theorems and lemmas, we show that a{Tt) is close to cr(T), which will be used in the 
proof of Theorem 13.11 

Theorem 4.6. Let Tt be the solution operator of the integral equation (1231), then 

a{Tt)C \J B{X^,2\\T-Tt\\H^^M)) (4-19) 

rieN 

Proof. Let tq = ||r — Tt\\H^(M)^ .4 = C\ UnGN B{Xn,2ro). For any z € A, using Lemma l4Tl we have 

\\{z-T)-^\\hi(m) < max I ^ 

riGN \z — Xn\ zrQ 

which implies that 

117" - Tt\\Hi{M) =ro< -T7 T-—TT| -■ 

Then using Theorem 14.11 we have z G p{Tt). 

Since z is arbitrary in A^ we get A C p{Tt). This means that 

a{Tt) = C\p{Tt) c C\.4 = IJ B(A„, 2||r - rt||^i(^)). 

nGN 


□ 

Now we can show Theorem 13.11 i.e., the convergence of the eigenproblem with the Neumann boundary . 

Proof. Let ri = ^ \\Tt,h - Ttllci + \\T-Uh^(mL -4 = C\ [U„eN 73(A„, r^) U S(0, . For any z G A, 

using Lemma 14.11 we have 

||(^-T,)-d|c^ < 

< 


(“is I" -^"1- -^‘11^0 


2|-A4| _ ||^_^|| V 

fk/A+l 11^ 7t||//l) 
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or 


as long as h small enough by using Theorem 14.41 
Both above two inequalies implies that 

\\Tt,h - TtW^, < 2||(^_T,)-i||cr 

Then using Lemma l4.11 we have z G p{Tt^h)- 

Since z is arbitrary in A, we get A C p(Tt^h)- This means that 

= CVm,^) C C\^ = y B(A„, n) U 5(0, (4.20) 

riGN 

Moreover, using Theorem 14.61 and the definition of ri, we have 

a(Ti) C U 5(A„,2ri). (4.21) 

kGN 


On the other hand, using Theorem 14.31 and 14.41 we know that there exist C > 0 independent on t and h, 
such that 


+ (4.22) 

For any fixed eigenvalue A„ £ cr(T), let 7 ^ = min |Aj — A|, j £ N and 7 = min 7 j. Using the structure 

AG<t(T)\{Aj} j<n 

of o-{T), we know that 7 > 0. Without loss of generality, we can assume t,h are small enough such that 
ri < 7 / 6 . 

Let Tj = {z £ C : |z — AjI = 7 / 8 }, Uj be the aera enclosed by Fj. Let 


<7tj- — cr{Tt) y Uj, o't,h,j — y Uj. 

Using the definition of F^ , we know for any j < n, Fj C p(T),p{Tt) and piTt^h)- 
In order to apply Theorem 14.21 we need to verify the conditions 


\\{T-Tt)Tt\\H^ < mm 


.er, ||( 2 -T)-i||ffU 


\\{Tt - Tt^h)Tt,h\\c^ < min 


\\{z-Tt)-^\cU 


(4.23) 

(4.24) 


Using Lemma 14.11 and the choice of Fj, we have 
kl > mincer, |z| 


— , 1 / ^ , 1 / —>(|Aj|- 7/3) min |z-A„| = (|Aj|- 7 / 3 ) 7 / 3 . 

zerjIK^: —T) ^\\hi max2grj IK^ “ ^) zGr^.nGN 


Then, using Theorem 14.31 and Lemma 14.51 condition 14.231 is true as long as t is small enough. 
Using Lemma 321 we have 


> 


mm^gr,- 


zerj 11(2; - Tt)-^\ci max^gr^ ||(^: - T't)"^llci 

(|Aj| -7/3)2#+2/4) 


> 


> 


2\M\ 

(|A,-|-7/3)¥^+^/")7 

12|7W| 


min | 2 ;-A„|-||r-rt||j^i 

z^Fj ,nGN 


(4.25) 


10 














To get the last inequality, we use the fact that \\T — Ti||j:^i < ri < 7/6 and min \z — A„| = 7/3. 

zGFj.neN 

Using Theorem 14.41 and Lemma lT5l we can choose h small enough such that condition 14.241 is satisfied. 
Then using Theorem 021 we have 

dim.{E{Xj,T)) = dim.{E{at,j,Tt)) = duji{E{at,h,j,Tt^h)), j = ,n. (4.26) 

Combining (I4.20[) . above equality would imply that 

|a7-A,|<c((‘'> + ^), ieN. (4.27) 

The convergence of eigenspace is also given by Theorem 14.21 For any x € £'(A„,T), ||a:||ci = 1, 

\\x-Eiat,n.Tt)x\\H^ < C\\T-Tt\\Hi\\x\\Hi<Ct^^^, (4.28) 

Ch 

\\E{at^^,Tt)x - E{at,h,n,TtMc^ < C m - + \\{Tt - Tt^h)TtA\c^) < (4-29) 

Finally, we have 

||x - E{at,h,u, Tt,h)x\\H^ < C . (4.30) 

□ 


5 Proof of Theorem 14.3 


We first introduce the local coordinate of the manifold M. According to Proposition 6.1 of HZ], M can be 
locally parametrized as follows. 

X = <I)(7) : fl C M'' -)■ A1 C (5.1) 

where 7 = ( 7 ^, • • • , 7 ^)* G and x = (x^, • • • , x'^Y £ M. 

Let di' = be the tangent vector along the direction 7 ® . Since A4 is a submanifold in with induced 
metric, 9^/ = (9^/$^, • • • , 9^/$'^) and the metric tensor 

gt'j' =< d^'.df >= 9i'$'9j'$'. 

Let < 7 * -1 denote the inverse of i.e.. 


gi'i' g 


I'f _ X. 

— Ui 


1, i = j, 
0, i Y 3- 


For any function / on Al, ^ f = g^ ^ dji fdi' denotes the gradient of /. For convenience, let V-’ / denote the 
A component of the gradient V/, i.e., 

= and 9,-/= 9,7$^^- (5.2) 


Then Theorem 14.31 is an easy corollary of following three theorems. 

Theorem 5.1. Assume A4 and dA4 are C°°. Let u(x) be the solution of the problem ()2.2I) and Ut{x.) be the 
solution of corresponding problem (12.3|) . If f £ C°°{M),g £ C°°{dM.) in both problems, then there exists 
constants C,Tq depending only on M. and dA4, so that for any t <Tq, 


Ltiu-ut) - / nj(y)77*(x,y)V*Vju(y)i?t(x,y)dTy 

J dM. 


V[Lt{u-Ut)- / nj(y)77i(x,y)V,Vju(y)At(x,y)dTy 
J dM 


L^{M) 


L^iM) 


where g{x,y) = ^^{x,y)di^{a) and a = ^ ^(x),^(x,y) = $ ^(x) — $ ^(y). 


< CY^^\\u\\h3{m)^ 

< C\\u\\h3{m)- 


(5.3) 

(5.4) 
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Theorem 5.2. Assume A4 and dM. are C°°. Let m(x) solves the integral equation 

LtU = r(x) — f (5-5) 

where r G H^{AA) and f = r(x)dx. Then, there exist a constant C > 0 ,Tq > 0 independent on t, 

such that 


Ikllifi(Al) < C (||r||i2(^) + <11 Vr||i2(^)) 

as long as t <Tq. 

Theorem 5.3. A ssume M. and dM are C°°. Let m(x) solves the integral equation 

Ltu= / bi{y)r],{yi,y)Rt{x,y)dTy -b 
J dJ^ 


where 77 is same as that in Theorem \5.1\ and 

1 


^=TT7t/ / ^*(y)^i(x,y)-Rt(x,y)dTydx. 


\M\ 


IM JdM 


Then, there exist constant C > 0 ,Tq > 0 independent on t, such that 


(5.6) 


(5.7) 


(5.8) 


as long as t <Tq. 

The proof of Theorem 15.21 can be found in m- In next two subsections, we will prove Theorem 15.31 and 
Theorem 15.11 sequentially. 


5.1 Proof of Theorem 15.31 

First, we need following four lemmas which have been proved in m- 

Lemma 5.1. For any function u G Lf{M), there exist a constant C > 0 independent on t and u, such that 

{u,Ltu)j^>C [ iVupdMx (5.9) 

JM 

where {f,g)_^ = /(x)5(x)d^x for any f,g G L 2 {M) and 

^^(x) = ^uiyWy, (5.10) 

wt(x) Jm \ J 

and u;t(x) = CtJj^R dfiy. 

Lemma 5.2. Assume M and dM are C°°. There exist a constant C > 0 independent on t so that for any 
function u G L 2 {M) with u = 0 and for any sufficient small t 

{Ltu,u)j^ > (5.11) 


Lemma 5.3. Assume both M and dM are C°° smooth. For sufficiently small t and any x G M, there are 
constants Wrain > 0,Wmax > 0 depending only on the geometry of M,dM, so that 


Wn 


< 


/ Rt{x,y)dpy<w 

n 

JM 


12 







Lemma 5.4. Assume both Ai and dAi are C°° smooth. For sujficiently small t and any x £ there 
is a constant C depending only on the geometry of Ai,dAi, so that for any integer s > 0 and i^t(x,y) = 
i?t(x,y),iit(x,y),^t(x,y), 

[ ^t(x,y)|x-y|M^y < [ ii:t(x,y)|x-y|®d/^y < 

JM JdM 

where Rt{x,y) = CtR and R{r) = /“ R{s)ds. 

Now, we are ready to prove Theorem 15.31 


Proof, of Theorem \5.3\ 

The key point is to show that 


'M 


t(x) 


'dM 


bii^WRt{yi,y)dTy - 6 d^; 


< CVt ||b(y)||iji(^)||'u||jYi(^) 


(5.12) 


Notice that 




IM JdM 

. Then it is enough to show that 


bziyHiA y)-Rt(x, y)drydx 


< CVi ||b(y)||i2(a^) < cVi ||b(y)||ffi(_M) 


u(x) 


>M 


'dM 


h{yW Rt{y^,y)ATy d/T; 


< Cy/t ||b(y)||/i-i(^)||M||^i(^) 


First, we have 


|2tV'^i?t(x,y) + Vi?t(x,y)| < C'|^|^i?t(x, y) 


(5.13) 


(5.14) 


where i?t(x,y) = CtR and i?(r) = /“i?(s)ds. 

This would gives us that 


[ u{x) f bi{y) (r]^Rt{x,y)+2t\7^Rt{x,y)) dTydfi, 
IM JdM ^ ' 


< C / |m(x)| / | 5 ,(y)||^| 2 i?t(x,y)dTyd/rx 


! M 


tdM 


<Ct |m(x)| / |6,(y)|i?t(x,y)dTydp, 


/M 


IdM 


<ct \Hy)\ / |'u(x)|i?t(x,y)drxdAiy 


IdM 


t M 


< Ct\\h{y)\\L2(^QM) ( 


\JdM \Jm 
2 


i?t(x,y)dp. 


'M 


i(x)|^i?t(x,y)dpx dr. 


1/2 


'dM 


< Ct\\h{y)\\Hi(M) / |u(x)| 

\Jm 

< C't^^'‘||b(y)||ji-i(^)||'u||i2(^) 

On the other hand, using Guass integral formula, we have 


i?t(x,y)dTy d/i; 


1/2 


(5.15) 


m(x) / 6i(y)Vii?t(x,y)dTyd/ix 


IM 


IdM 


IdM 


bi{y) / M(x)Vii?t(x,y)d/ixdry 


IM 


bi{y) ni(x)u(x)i?4(x,y)dTxdTy- 


divx[&i(y)u(x)]i?t(x,y)d/ixdry. (5.16) 


IdM 


IdM 


IdM JM 
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For the first term, we have 


/ bi{y) / nj(x)'u(x)i?t(x,y)dTxdry 
' dM. J dM. 


< / l^*(y)l / |u(x)|i?t(x,y)dTxdTy 

J dM. J dM 

< C'||b(y)||i 2 (a_M) ( f ( f |u(x)|Et(x,y)drx) dry ) 

V JdM dM / } 


'dM 


< C'l|b(y)||i/i(Ai) / / i?t(x,y)dT, 

\*/ dM \JdM 

< ||b(y)||ffi(7K) f / (f ^t(x,y)dry^ |M(x)pdrx^ 

dM \J dM J J 


|M(x)pi?t(x,y)drx dr. 


1/2 


1/2 


< Ct ||b(y)||jYi(^)||u||i 2 (a^) 

< ||b(y)||_H-i(x)||u||fl-i(;v(). 

We can also bound the second term of (15.161) 


(5.17) 


IdM JM 
< C'l|b(y)||L2(a^) 


divx [ 6 i (y )m (x)] (x, y) d/ix dry 

I V'u(x) I .Rt (x, y) d/ix dry 


IdM JM 


< C'l|b(y)||ffi(7W) ( ^t(x,y)dMx ) dry 

■ 

< C'l|b(y)||ffi(Ai) ( / |Vu(x)|M / i?t(x,y)dry ) d/i: 

VJAt VdOAt / 

< ||b(y)||/i-i(Ai)||u||fl-i(A^)- 


1/2 


(5.18) 


Then, the inequality (15.131) is obtained by (15.151) . (I5.16L (15.171) and (15.181) . 
Now, using Lemma 15.21 we have 


Il 2 (A 4 ) ^ C' {u, Ltu) — C 


/ ■u(x) / 6*(y)77'^i?t(x,y)dryd/ix 

'M JdM 


< CVt ||b(y)||/i-i(^)||M||/i-i(^). (5.19) 


Denote p(x) = 6 i(y)? 7 '’i?t(x, y)dry, then direct calculation would give us that 

lb(x)|U2(;^,) < C'b/'^||b(y)||^^i(^), 
l|Vp(x)||i 2 (;^,) < Ct-^^^\\h{y)\\Hi(M)- 

The integral equation LtU = p gives that 


u(x) = u(x) + 


t 


Wt{-x.) 


p(x) 


where 


K^) = —^[ -Rt(x,y)M(y)d/iy, wt{x)=[ i?t(x,y)d/iy 
Wtw Jm Jm 


(5.20) 

(5.21) 


(5.22) 


(5.23) 
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Then by Lemma l5.11 we have 

W^uWl.^M) < 2||Vn||i.(^)+2i2 


P(x) 

Wt(x) 


L^{M) 

2|ir7^ll2 


< C {u, Ltu) + Ct\\p\\j^ 2 (^_f^'j + Ct ||Vp||^ 2 (^) 

< C\/t ||b(y)||/i-i(^)||u||//i(^) + C't||p|1^2(;v^) + Vp|1^2(;v^) 

< C'l|b(y)||_f/i(Ai) (v^||'«||j/i(A^) ■ 

Using (15.191) and (|5.24l) . we have 

\W\\‘h^{m) - C'||b(y)||//i(;v^) (yt\\u\\H^{M) 

which proves the theorem. 


(5.24) 

(5.25) 

□ 


5.2 Proof of Theorem 15.11 

Proof. Let r(x) = LfU — L^ut^ then we have 

r(x) = -j [ i?t(x,y)(M(x)-u(y))dAiy+2 [ Rt{x,y)g{y)dTy - ( Rt{x,y)f{y)dpy 
£ Jm JdM J M 


i?t(x, y)(M(x) - u{y))dny + / i?t(x, y)AMu{y)dny 


IM 


IM 


+ 7 / (x-y) • Vu(y)i?t(x,y)dd^j 
£ Jm 


t 


i?t(x, y)(-«(x) - u(y) - (x - y) • V'u(y))d/ry + / i?t(x, y)AMu{y)dpy 


(5.26) 


<M 


IM 


Here we use that fact that u is the solution of the Poisson equation with Newmann boundary condition (12.21) , 
such that 


IdM 


f du 

■Rt(x,y)3(y)dTy = / i?t(x,y)^(y)dry 

JdM £££1 


and 


Denote 


/ i?i(x,y)A7ciM(y)d/ry + / Vyi?t(x, y)Vu(y)dAiy 
IM Jm 

/ Rt{^,y)^Mu{y)A^ly + - (x-y) • V'u(y)i?t(x,y)dd/ry, 


/ Rt{^,y)f{y)dny = Rt{x,y)AMu{y)dpy. 
IM Jm 


S; = {y e 7W : |x - y| < r} , = {y e : |x - y|^ < 32t} 


(5.27) 


where r > 0 is a positive number which is small enough. 

Since the manifold is compact, there exists x^sAf, i = 1, - ■ ■ ,N such that 


N 


Me 


(5.28) 


By Proposition 6.1 in we have there exist a parametrization : fii C —>■ Ui C M,; i = 1, • • • ,JV, 

such that 
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1. C Ui and is convex. 

2. $ G C^{n)- 

3. For any points x,y G rj, ^ |x - y| < ||$i(x) - $*(y)|| < 2 |x- y|. 

Moreover, we denote $(/3) = x, $(a) =y,^ = P — a,r] = and 

ri(x) = -i ^ i?t(x, y) ^m(x) --«(y) - (x-y) • Vu(y) - i77*77'’(V*V'^u(y))^ d/Ty 

?'2(x) = 4 / i?t(x,y)77*7?'^(V*V%(y))dAiy - [ u{y)V^Rt{x, y)dny 

^ 3 (x) = [ r]\V^V^u{y)V^Rt{x, y)dny + [ div {r]^u{y)) Rt{x, y)dfiy 

Jm Jm 

r 4 (x) = [ div (77*(V*V'^u(y)) .Rt(x,y)d//y + / RtiT<i,y)AMu(y)dny. 

Jm Jm 


then 


r(x) = ri(x) - ri(x) - 7-3 (x) + r4(x). 

Next, we will prove the theorem by estimating above four terms one by one. 
First, let us consider ri. To simplify the notation, let 

1 


Then, we have 


d(x, y) = u(x) - u{y) - (x - y) • V'u(y) - -77*77^ (V*V^M(y)). 


/ \ri{x)fdn^ = / Rt{x,y)d{x,y)dfiy d^^ 

J A 4 J Ai J Ai 

i?t(x,y)d7iy^ (^J i?t(x,y)|d(x,y)pdAiy^ d/r^ 


M \JM 


< C I I i?t(x,y)|d(x,y)pd 77 yd 7 ix 

■Rt(x,y)|d(x,y)pd7iyd/ix 


JM JM 

j 

Jm Jb 

N 


! JBC 


i?t(x,y)|d(x,y)pd7iyd/rx 


= ^E 

= f f f d/T; 


(5.29) 


(5.30) 


(5.31) 
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Using the fact that is convex and the Newton-Leibniz formula, we can get 


d(x,y) = u(x) - u(y) - (x-y) • VM(y) - 


^ I I I + S 3 S 2 S 1 C))) ds 3 ds 2 dsi 

= CC'C" f f f sls2d,<^^ {a + S3Si^)d,>'di><^^' {a + S3S2SiOV^' V^u{<^{a + S3S2Si^))ds3ds2dsi 

Jo Jo Jo 

+CC C [ [ [ sldi»di^^{a +S3Si^)di>^^ {a +S3S2Si^)V^ V^u{^{a +S3S2Si^))ds3ds2dsi 

Jo Jo Jo 

+CC C f f f sls2di^^{a +S3S2Si^)di>^^ {a +S3S2Si^)di»^^ (a + S3S2Si^) 

Jo Jo Jo 

V-^ V-^ V'^lt(<I)(Q! + S 3 S 2 Sl^))dS 3 dS 2 dsi 
Using this equality and $ G C^(fl), it is easy to show that 

f f [ f i ^t(x,y) |-D^’^u(d>i(a + S 3 S 2 Si^))|^d/ixd/Xyds 3 ds 2 dsi 
Jo Jo Jo Jb^. Jm^. 

max [ [ i?t(x,y) + s^))|^d/rxd/Xy (5.32) 

0 <s<l Jb;. dMJ, 

where 

n n 

\D^'^u{x.)\^ = |V^"v^''v^u(x)| 2 + ^ |V^'v^u(x)| 2 . 

Let Zi = $i(a + s^), 0 < s < 1, then for any y G B^. and x G A^y, 

|zi - y| < 2 s|^| < 4s|x - y| < 8sVi, |zj - x*] < |zj - y| + |y - x^l < r + 8 sv^. 

We can assume that t is small enough such that 8^/i < r, then we have 

z* e 

After changing of variable, we obtain 

/ / i?t(x,y) |D^’^M(4>j(a + S 3 S 2 Si^))|^d/ixd/iy 

dsr Jmi 


< 


< 


< — 


^0 Jb^ Jb^-^ s 

C 

So JB^. S 


] |D^’^'u(zi)|"^ d^z-d^j 


128s2t 


LuflELZll. 

[ 128 sH 


^ dfiy J |£)^’^u(zi)|^d/rz 


< C f |D^’^u(x)|^ d/ix 

■’Bi: 

This estimate would give us that 

lkl(x)||L2(;K) < Ct^/^u\\H3(M) 


(5.33) 


(5.34) 
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Now, we turn to estimate the gradient of ri. 

[ |Vxri(x)pd^x < C f [ Vxi?t(x,y)d(x,y)d/ry 

Jm Jm Jm 


d/^x C 


fM 


i?t(x,y)Vxd(x,y)dAij 


'M 


d/Tx- 


Using the same techniques in the calculation of ||ri(x)||i 2 (^), we can get that the first term of right hand 
side can bounded as follows 


Jm 


< M 


Vxi?t(x,y)d(x,y)d/rj 


d^x < C\\u\\'jj3(^j^y 


(5.35) 


The estimation of second term is a little involved. First, we have 


JM 


■Rt(x,y)Vxd(x,y)d/rj 


'M 


d/Tx 


< 


i?t(x,y)dAiy^ (^J i?t(x,y)IVxd(x, y)pd/Ty^ d/Tx 

^i(x, y)|Vxd(x, y)pd/ryd/Xx 


JM \JM 

< c [ [ ^ "^'2 

Jm Jm 

< CV [ [ i?t(x,y)|Vxd(x,y)pd^yd/ix 

Jm Jb-^ 

N ^ /• 

^t(x, y)| Vxd(x, y)|^d^yd/i: 


Jb^’: J b^ 


= ^E 

2 = 1 ' 

= C'E f if ^t(^>y)l^xd(x,y)pdAix ] d/Ty. 
z=i Jb-^ \Jm^^ J 


(5.36) 


Using Newton-Leibniz formula, we have 

d(x,y) = CC J J Si {a + {a + S 2 Si^)V^ + S 2 Si^))j ds 2 dsi 

-ff ^ i) '51 ds2dsi 

Then the gradient of d(x, y) has following representation, 

Vxd(x,y) = ^ si(5,$i(a + si05*^$i'(a + S2Si0Vi'vin($(a + S 2 Sie)))ds 2 dsi 


(5.37) 


(^CC ^ J J J + sssiO^i'd)! (a + S3S2SiC)Vi V'iM(<i)(Q; + S 3 S 2 S 1 O)) ds3ds2dsi 

(5.38) 


+V, 


/O ^0 ^0 

= di(x,y) +d 2 (x,y). 


For di, we have 


< 


Ib (/ * 

Ct'^ [ [ [ if i?t(x,y)|i:)^’^w($(a + S 2 SiO)pd/ix I d/iyds 2 dsi 

Jo Jo Jb;^. \J-My J 

^ i?t(x,y)|i:)2’3u($(a + sO)Pd/rx^ d/Ty 


< Ct^ max 

0<s< 


(5.39) 
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which means that 


[ if -Rt(x,y)|di(x,y)pdAix ) d^iy < C / \D'^’^u{x)\'^dn^ 

JB^, J JBi': 


(5.40) 


For 6 , 2 , we have 
rf2(x,y) 

= Vx (^CC ) J (« + S3S2Si^)V'^ V%($(a + S3S2S1C))) ds3ds2dsi 


= V: 

+ V 
+ V 


(«■')«■" J 

(CF) F 

(«■«■■) C" 


sls2^^^^{a + S3sy)de>di>^^ {a + S3S2SiC)V^ V'^u( 4 >(a + S3S2Si^))ds3ds2dsi 


'[ 0 . 1 ] 3 


{a + S3Si^)di>^^ {a + S3S2sy)W^ V''zi($(a + S3S2Si^))ds3ds2dsi 


'[ 0 . 1 ] 3 


s?S 29 jd>''(a + S 2 SiC) 9 j' 4 >^ (a + S 3 S 2 SiC) 9 j" 4 >'^ (a + S3S2S1C) 


V-’ V-^ V-^w(<i)(Q! + S3S2Si^))ds3dS2dsl 


This formula tells us that 


[ if ^t(x,y)|d2(x,y)pd^x I d/iy 
Jb-^ \Jm^^ j 



Using the same arguments as that in the calculation of Ilri||i 2 (^), we have 


f if -Rt(x,y)|d 2 (x,y)pd/ix I d/iy < C / |D^M(x)pd/rx 

JBJ. / JBl'-, 


Combining (I5.40p and (15.421) . we have 


I^^’i(x)||l2(ai) < C'II'w||b3(ai) 


(5.41) 


(5.42) 


(5.43) 


The estimates of r 2 , r 3 and are similar as those in our previous paper m- In order to make this proof 
self-consistent, we also give a complete proof of this part. 

For r 2 , first, notice that 

V^Rt{x,y) = ^am'4>'’(a)5™'"'9„'4>*(a)(a;* - y*)i?t(x,y), 

^i?t(x,y) = ia™,4>^(a)5™'"'a„-$*(a)ea,,$»i?t(x,y). 

Then, we have 


V'^i?t(x,y) - ■^i?t(x,y) 

= - yn - Vi?t(x,y) 

^ - C 9, i?t(x,y) 




sdj>di>d>^{a + rs^)dTds ) iit(x, y) 


'0 Jo 
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Thus, we get 


Rtix,y) - ^iit(x,y) 




Vx ( V'^i?t(x,y) - —i?t(x,y) 




Then, we have following bound for r2, 

[ k2(x)pd^x < ct f f f Rt{x,y)\D'^u{y)\d^j,y\ d^x 
jm Jm \Jm J 

< Ct f ( f i?t(x,y)d/iy) [ Rt{x,y)\D'^u{y)fdfiydn^ 

Jm \Jm / Jm 

< Ctmax^y i?t(x,y)d/ix^ J \D'^u{y)\'^d^iy 

Similarly, we have 

f |Vr 2 (x)pd/ix < Ct [ ( f Vxi?t(x,y)d/ry^ [ Vxi?t(x,y)|-D^M(y)pd/ryd/rx 
Jm Jm \Jm J Jm 

< CVtmax^y Vxi?t(x,y)d/rx^ max ^y Vxi?t(x,y)d/ry^ J \D'^u{y)\'^dfi 


(5.44) 




y 

(5.45) 


ra is relatively easy to estimate by using the well known Gauss formula. 

r3(x) = f n^r]\V'-y^u)Rt{yi,y)dTy (5.46) 

J d^A 

Now, we turn to bound the last term r^. Notice that 

VJ(V^u) = dv {du'U)) 

= {dv{d^.^^)) g^'^\dn'u) + {du'^^)g^'^' (ff”^'"'(5„m)) 

= -^{dm'y/g)g^'^'{dn'U)+dm' (5™'"'(5„m)) 

= (V5g™'”'(a„-u)) 

= I^MU. (5.47) 


Here we use the fact that 


= {dm'{dk'<^^))g^'^' {dv^=) 

= \g^'^'dm'igk'v) 

= -^{.dm’y/g) (5.48) 

Va 
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Moreover, we have 


= -Z^' {dy ) {di> $*) (V* u) 

= (v%) 

= {dy^^)d,, {V^u) 

= -V^ (V%) . (5.49) 

where the first equalities are due to that di'^^ = —5\,. Then we have 


div (77*(V®V^M(y))) + A^it 

= -^d,, (ZZ/^'(a/4>^)e'(5i$*)(v*v%(y))) -/^■'(a,v$^)(a,,e')(az$*)(v*v%) 

Here we use the equalities (15.471) . (I5.49L rf = and the definition of div, 

1 


divX = 


'JWi 


dA^W'^' dy<^^X^). 


(5.50) 


(5.51) 


where A is a smooth tangent vector field on Xi and is its representation in embedding 

coordinates. 

Hence, 


r4(x) 




(v^/^'(9,.4>^)(9i$*)(V*V%(y)))i?t(x,y)d/iy 


Then it is easy to get that 


lk4(x)||L2(;^,) < Ct^/'^\\u\\H^(^M), 

||Vr4(x)||^2(;^) < C\\u\\h3{m)- 

The proof is complete by combining (15.34^ . (15.431) . (15.44|) . (15.45^ . (15.461) . (153^ . (15.531) . 


(5.52) 

(5.53) 

□ 


6 Proof of Theorem 14.4 

First, we need the following theorem which have been proved in m- 

Theorem 6.1. Assume M and dM are C°° and the input data (H, V) is an h-integral approximation of 
M.. Let f € C(Ai) in both problems, then there exists constants C,To depending only on M. and dM. so that 

Ch 

\\Tt,hf-T,f\\H^M) < ^ll/lloo, (6.1) 

as long as t < Tq and < Tq. 

Then the main idea to prove Theorem 14.41 is to lift the covergence from to by using the fact that 
Ttu and Ttyu have higher order regularity for any u € C(M). The details are given as following. 

Proof, of Theorem 
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First, for any u £ C^{Ai), we know that Ttu and Tt^hU have following representations 


Ttu = f i?t(x,y)TtM(y)dy H-^ [ i?(x, y)'u(y)dy, 

wt(x) wt(x) 

Tt,hU = -;^y^i?t(x,Xi)Tt /jM(xi)t/i H- ^-^y^ R{x.,:x.i)u{-Xi)Vi. 

... ^ wt^hm ^ 


Wt,hi^) 


Denote 


T/it = 


T?u = 


Wt,/t(x) Jj^ 


Wt,h{^) Jm 


i?t(x,y)rtM(y)dy H- [ i?(x, y)M(y)dy, 

Wt^hy^) Jji4 

i?t(x,y)rt,?,M(y)dy + / ^(x, y)M(y)dy, 


wt,h(x) Jm 


Direct calculation would give that 


||■u;^(x) - ■u;t,?>(x)||oo < ||Vu;t(x) - Vwt,ft(x)||oo < ^ 


and 


IM 


Rt{^,y)Ttu{y)dy 


Vx [ i?t(x,y)rtM(y)dy 
Jm 


< C't-('=+2)/4||Ttw|U2 < Ct-('=+2)/4||w||ci 


and 


'M 


i?(x, y)'u(y)dy 


<C||n|U, 


Vx / i?(x,y)u(y)dy 
Jm 




Using above inequalites, we have 


\TtU — TIu\ < 


Ch 


which proves that 


Secondly, we have 


i “I - i(fc+2)/4l 
Ch 

\V{T.u-Tlu)\ < 




T}u-T?u\ = 


wt,h{:!i-) Jm 


Rt (x, y) (Ttu(y) - Tt^hu[y)) dy 
Ch 




and 


\V{Tju-Tju)\ = 


wt,h{yj Jm 


-Rt(x, y) {Ttu{y) - Tt,hu{y)) dy 
Ch 




( 6 . 2 ) 

(6.3) 

(6.4) 

(6.5) 

( 6 . 6 ) 

(6.7) 


( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 


(6.13) 


(6.14) 
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This implies that 


Tj M Tj U pi < ^fc/4_,_2 1'“ loo 

(6.15) 

Using Theorem 14.51 we have 


\Tt,hu\ < Ct-'=/4||u|U, 

(6.16) 

m,hu\ < Ct-('=+2)/4||u|U. 

(6.17) 

Thus, 


\Tt^hU Tj u| < ^(■jij^2)/4 ll'^ll°°’ 

(6.18) 

|V(rp,u-T2u)| < ^^llulloo. 

(6.19) 

which also reads 


||Ft,/iW Ft — ^fc/4+l 11^11°° 

(6.20) 

The proof is complete bv combining (|6.12p. (|6.15p and (|6.20p. 

□ 


7 Proof of Theorem 14.5 

In this section, we prove the bounds on the solution operators Tt and Tt^h- First, we need following theorem 
whose proof can be found in m- 

Theorem 7.1. If both the submanifolds M and dAi are C°° smooth and the input data (P, S', V, A) is an 
h-integral approximation of then there exist constants C > 0, Co > 0 independent on t so that for 

any u = (ui, • • • , M|p|)* with X]l=i = 0 and for any sufficient small t and 

(u, Tu)^ > C(1 - (u, u)v , 

where (u,v).^ = Y}^[utViVi for any u= {ui,--- ,M|p|)‘,v = (ui,--- ,U|p|)‘. 

Now, we have enough materials to prove Theorem 14.51 
Proof, of Theorem \4-.5\ 

Let u = Tt{f). From Lemma [5.21 there exists a constant C independent of t so that 

< c'II'w||l2(ai)||/||l2(a^), 

which means ||w||p 2 (^) < C||/||p 2 (^). At the same time, we can write 

w(x) = —^ / Pt(x,y)M(y)d/iy-Pt(x,y)/(y)d^y. 

Wt(^) Jm Wt^xj jj^ 

Since the kernel function, R, is chosen to be C^, it is also belong to C^ and 

|m(x)P < c( [ Rt{x,y)dp.y] ( f Pt(x,y)-«^(y)dAiy 

\JM / \Jm 

+Cf (^J^Rt{x, y)dny^ ^^Pt(x,y)/2(y)d/ry^ 

< CCt (||u||i2(^) +t^||/||i2(^)) < Ca||/||i.(^) 
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Similarly, we have 


||Vm||oo < ^(I|m||oo +t||/||oo) < 

Now, we turn to bound Tt^h- We can write 


^^1/2 

^ll/lloo. 




1 

Wt,h{^) 


M^,P3)UjVj 

Pj^p 


t 

Wt,/i(x) 


i?t(x,Pj)/(pj)Vj-, 

Pj&p 


where u = (ui, • • • , U|p|)* with X]l=i = 0 solves the problem (12.51) . 

From Theorem O we can easily get that — C'||/||oo for some constant C independent 

on t. It is obvious that for some constant C independent on t, 

/|P| 

\TtAf)i^)\ < ccY I I +^*ll/lloo < ccYwfiA, 

and 

1 /2 

^^ 1/2 / | P | \ ^^ 1/2 

□ 


8 Conclusion and Future Work 


In this paper, we proved that the convergence of the point integral method for the spectra of the Laplace- 
Beltrami operator on point cloud. And the rate of convergence is also obtained. This work builds a solid 
mathematical foundation for many Laplace spectra based algorithm. 

In many applications, the sample points of the manifold are draw according to some probability distribu¬ 
tion. Then one interesting problem is to study the performance of the point integral method on the random 
samplesas the number of sample points tends to infinity. Based on the results reported in this paper We can 
show that with overwhelming probability, the spectra given by the point integral method converge to the 
spectra of following eigen problem 

I X e 

1 §^(x)=0, xedM, 


where p is the probability distribution. The rate of convergence can also be obtained. This result will be 
reported in our subsequent paper. 
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